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Abstract 

Qs ■ We evaluate the infinite volume, continuum limit of glueball masses in the 

valence (quenched) approximation to lattice QCD. For the lightest scalar and 
-J , tensor states we obtain masses of 1648 it 58 MeV and 2267 it 104 MeV, re- 

Q I spectively. 
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I. INTRODUCTION 

In recent articles we described calculations of the infinite volume, continuum limit of 
scalar and tensor glueball masses in the valence (quenched) approximation to lattice QCD 
IHJ^ . For a single value of lattice spacing and lattice volume, we reported also a calculation 
of the the decay coupling constants of the lightest scalar glueball to pairs of pseudoscalar 
mesons. The mass and decay calculations combined support the identification of /o(1710) 
as primarily composed of the lightest scalar glueball 0,§]. Evaluation of the mass of the 
lightest scalar quarkonium states and of quarkonium-glueball mixing amplitudes f^ then 
yield a glueball component for /o(1710) of 73.8 ± 9.5 %. In the present article, we describe 
the glueball mass data of Ref. |Q in greater detail along with an improved evaluation of the 
mass predictions which follow from these data. For the scalar and tensor glueball masses we 
obtain 1648 ± 58 MeV and 2267 ± 104 MeV, respectively. 

The valence approximation, on which our results depend, may be viewed as replacing 
the momentum dependent color dielectric constant arising from quark-antiquark vacuum 
polarization with its zero-momentum limit and, for flavor singlet mesons, shutting off 
transitions between valence quark-antiquark pairs and gluons. The valence approximation 
is expected to be fairly reliable for low lying flavor nonsinglet hadron masses, which are 
determined largely by the low momentum behavior of the chromoelectric field. This expec- 
tation is supported by recent valence approximation calculations [^|[ of the masses of the 
lowest flavor multiplets of spin 1/2 and 3/2 baryons and pseudoscalar and vector mesons. 
The predicted masses are all within about 10% of experiment. For the lowest valence ap- 
proximation glueball masses, the error arising from the valence approximation's omission 
of the momentum dependence of quark-antiquark vacuum polarization we thus also expect 
to be 10% or less. Refs. 0,0 show this error should tend to lower valence approximation 
masses below those of full QCD. For flavor singlet conflgurations whose quantum numbers, 
if realized as quarkonium, require nonzero orbital angular momentum, it is shown in Ref. 
that the additional error arising from the valence approximation's suppression of transitions 
between valence quark-antiquark pairs and gluons is likely to introduce an additional error 
of the order of 5% or less. For the lowest scalar glueball this error is examined in detail in 
Ref. and found to shift the valence approximation mass by about 5% below its value in full 
QCD. It is perhaps useful to mention that, for glueball masses, the valence approximation 
simply amounts to a reinterpretation of the predictions of pure gauge theory. 

In Section || we deflne a family of operators used to construct glueball propagators. 
In Section |ITI| we describe the set of lattices on which propagators were evaluated and the 
algorithms we used to generate gauge conflgurations and estimate error bars. In Sections |^ 
and |V| we present our results for scalar and tensor glueball propagators, respectively, and 



masses extracted from these propagators. In Section we estimate the difference between 



the scalar and tensor masses we obtain in flnite volumes and the corresponding inflnite 



volume limits. In Section VII we extrapolate scalar and tensor masses to their continuum 



limits. In Section Vlll we compare our calculations with work by other groups [n],rT|. For 



combined world average valence approximation scalar and tensor glueball masses we obtain 
1656 ± 47 MeV and 2302 ± 62 MeV, respectively. 



II. SMEARED OPERATORS 

We evaluated glueball propagators using operators built out of smeared link variables. 
Glueball operators built from link variables with an optimal choice of smearing couple more 
weakly to excited glueball states than do corresponding operators built from unsmeared 
links. As a consequence, the plateau in effective mass plots for optimally smeared operators 
begins at a smaller time separation between source and sink operators, extends over a larger 
number of time units, and yields a fitted mass with smaller statistical noise than would be 
obtained from operators made from unsmeared link variables. Examples of the improvements 
which we obtained by a choice of smeared operators with be given in Section |rV| and |V[ 

Initially, we constructed smeared operators from gauge links fixed to Coulomb gauge. 
This method gave adequate results for the largest values of lattice spacing we considered. 
As the lattice spacing was made smaller, however, we found that the computer time required 
to gauge fix a large enough ensemble of configurations to obtain useful results became 
unacceptably large. We then switched to a gauge invariant smearing method. For the lattice 
sizes used in our extrapolation to the continuum limit, the gauge invariant mass results had 
statistical uncertainties typically a factor of three smaller than our earlier Coulomb gauge 
results. In the remainder of the present article we discuss only the gauge invariant results. 
A summary of our Coulomb gauge mass calculations is given in Ref . |]I|] . 

A family of gauge invariant smeared operators we construct following the adaptation in 
Ref. |T^ of the smearing method of Ref. []T3[. A related method of gauge invariant smearing 



is proposed in Refs. |T^. For ra > 0, e > 0, we define iteratively a sequence of smeared, 
space-direction link variable Uf'^{x), with ?7°'^(x) given by the unsmeared link variable Ui{x). 



Let Ui (x) be 
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eY.[uri^-j)]^ur{x-j)ur{x+z-j), (2.1) 
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where the sum is over the two space directions j orthogonal to direction i. The projection 
of ul^ (x) into SU(3) defines the new smeared link variable U^ (x). 
To find [/j (x) we maximize over SU(3) the target function 

ReTr{ut'''^\x)[u^^'^\x)]^}. (2.2) 

The required maximum is constructed by repeatedly applying an algorithm related to the 
Cabbibo-Marinari-Okawa Monte Carlo method. We begin with f/j (x) chosen to be 
1. We then multiply f/j (x) by a matrix in the SU(2) subgroup acting only on gauge 
index values 1 and 2 chosen to maximize the target function over this subgroup. This 
multiplication and maximization step is repeated for the SU(2) subgroup acting only on 
index values 2 and 3, then for the subgroup acting only on index values 1 and 3. The entire 
three step process is then repeated five times. Five repetitions we found sufficient to produce 
a ?7j (x) satisfactorily close to the true maximum of the target function in Eq. ( p. 21) . 
Iteratively maximizing the target function over SU(2) subgroups turns out to be much easier 
to program than a direct maximization over all of SU(3). The additional computer time 



required for this iterative maximization, on the other hand, was a neghgible fraction of the 
total time required for our calculation. 

From the Ul"'{x) we construct WJ^^'^^x) by taking the trace of the product of U^^ix) 
around the boundary of an s x s square of links beginning in the k direction. The sum 
of W^f^{x) over all sites with a fixed time value t gives the zero- momentum loop variable 

For each triple (n, e, s), a field coupling the vacuum only to zero-momentum scalars states 
is 

S"^^{t)=Y^ReWi;^%t), (2.3) 

where the sums are over space directions i and j. A possible choice of the two independent 
operators coupling the vacuum only to zero-momentum tensor states is 

T[^"(t) = 2ReWl'^'%t) - ReW^^'it) - ReW^^''{t) 

T^'^it) = V3ReW^:^''{t) - V3ReW^^''{t). (2.4) 

The optimal choice of (n, e, s) for each operator and lattice spacing will be considered in the 
next section. 



III. LATTICES, MONTE CARLO ALGORITHM AND ERROR BARS 

The set of lattices on which we evaluated scalar and tensor glueball propagators is listed 
in Table |. 

On each lattice, an ensembles of gauge configurations was generated by a combination of 
the Cabbibo-Marinari-Okawa algorithm and the overrelaxed method of Ref . |]15| . To update 



a gauge hnk we first performed a microcanonical update in the SU(2) subgroup acting on 
gauge indices 1 and 2. This was then repeated for the SU(2) subgroup acting on indices 
2 and 3, and the subgroup acting on indices 1 and 3. These three update steps were then 
repeated on each link of the lattice. After four lattice sweeps each consisting of the three 
microcanonical steps on each link, we carried out one Cabibbo-Marinari-Okawa sweep of the 
full lattice. 

At least 10,000 sweeps were used in each case to generate an initial equilibrium con- 
figuration. The number of sweeps skipped between each configuration used to calculate 
propagators and the total number of configurations in each ensemble are listed in the third 
and fourth columns, respectively, of Table |l|. Although the number of sweeps skipped in 
each case was not sufficient to permit successive configurations to be treated as statistically 
independent, we found successive configurations to be sufficiently independent to justify the 
cost of evaluating glueball operators. 

For the propagators, effective masses and fitted masses to be discussed in Sections |^ 
and|V], we determined statistical uncertainties by the bootstrap method [|T6|. The bootstrap 



algorithm can be applied directly, however, only to determine the uncertainties in quantities 
obtained from an ensemble whose individual members are statistically independent. We 
therefore partitioned each ensemble of correlated gauge configurations into successive dis- 
joint bins with a fixed bin size. Bootstrap ensembles were then formed by randomly choosing 



a number of entire bins equal to the number of bins in the original partitioned ensemble. 
For bins sufficiently large, propagator averages found on distinct bins will be nearly inde- 
pendent. It follows that for large enough bins, the binned bootstrap estimate of errors will 
be reliable. It is not hard to show that once bins are made large enough to produce nearly 
independent bin averages, further increases in bin size will leave bootstrap error estimates 
nearly unchanged. The only variation in errors as the bin size is increased further will come 
from statistical fluctuations in the error estimates themselves. To determine the required bin 
size for a particular error estimate to be reliable we applied the bootstrap method repeatedly 
with progressively larger bin sizes until the estimated error became nearly independent of 
bin size. The final bin size we adopted for each lattice, chosen to be large enough for all of 
the error estimates done on that lattice, is given in fifth column of Table |l[ 

IV. SCALAR PROPAGATORS AND MASSES 
From the scalar operator of Eq. |2.3| , a propagator for scalars is defined to be 

Pn^i - ^2) =< 5""(ti)5""(t2) > - < 5""(ti) >< 5""(t2) > . (4.1) 

To reduce statistical noise, Pg^'^{ti — 12) is then averaged over refiections and time direction 
displacements of ti and ^2- 

The collection of values of smearing iterations n, smearing parameter e, and loop size s 
for which propagators were evaluated for each lattice are given in Table |I|. At (3 of 5.70, 
and at j3 of 5.93 on the lattice 16^ x 24, we ran with relatively larger ranges of parameters to 
try to find values which coupled efficiently to the hghtest scalar glueball. For other lattices, 
the parameter range was then narrowed to choices which, in physical units, were about the 
same as the value range which gave best results at 13 of 5.93. 

From the existence of a self-adjoint, positive, bounded transfer matrix for lattice QCD, 
it follows that a spectral resolution can be constructed for Pg'^'^{t), 

Pnt) = E Z,{exp{-E,t) + exp[-E,{L - t)]}, 

i 

Zi = \< i\S''''{0)\vacuum > |^ (4.2) 

where the sum is over all zero- momentum, scalar states < i\, Ei is the energy of < z|, and L 
is the lattice period in the time direction. For large values of t and L, the sum in Eq. ( |4.2|) 
approaches the asymptotic form 

P^"(t) -^ Z{exp{-mt) + exp[-m{L - t)]} (4.3) 

where m is the smallest Ei and thus the mass of the lightest scalar glueball and Z is the 
corresponding Zj. Fitting Pg'^^{t) to the asymptotic form in Eq. ( |4.3| ) at t and t + 1 gives 
the scalar effective mass m(t), which at large t approaches m. 

To extract values of m from our data sets, we began by examining effective mass graphs 
to find combinations of n and s for which m{t) shows a plateau at t values for which we 
have data, and to determine which of these combination of n and s has the best plateaus. 

Among the data sets used in our final extrapolation of the scalar mass to zero lattice 
spacing, we included the largest range of values of n and s for the lattice 16'^ x 24 with /3 



of 5.93. Scalar effective masses obtained for this case with n of 5 and s of 3 — 7 are shown 
in Figures |l| - ^, respectively. As the loop size s is increased, initially the effective mass 
graphs become flatter, as shown, for example, by a decrease in the difference between m{0) 
and m{2). It follows that the relative coupling of the corresponding operators to the lightest 
scalar glueball increases with s. Beyond s of 5, however, this trend reverse. Thus, as might 
be expected, the relative coupling to the lightest state becomes weaker again when the loop 
is made too large. For s of 7 the effective mass graph shows no sign of becoming flat even at 
the largest t for which we have statistically significant data. For n of 5, the best coupling to 
the lightest state appears to occur with s of 4 or 5. For s fixed at 4, Figure ^ and Figures |^ 
- ^ show the variation in the effective mass graph as n runs from 5 to 8, respectively. The 
difference between m(0) and m{2) is least at n of 6 and then grows again as n is raised 
toward 8. 

In Figures [l| - H the statistical uncertainty in effective masses grows as t is made larger 
and tends to grow also if n or s is increased. Both of these phenomena are explained by the 



discussion in Ref. [[T^l of the statistical uncertainty in propagators. 

Figures |^ - show scalar effective masses for the each of the values of lattice size and 
(3 listed in Table |. The parameters n and s for the data in Figures |^ - |13| are chosen, for 
each lattice and jS, from among the set which couple best to the lightest scalar. 

For each combination of lattice size and (3, we determined a final value of the scalar 
mass from the collection of propagators for which the effective mass graph showed at least 
some evidence of a plateau at large t. For several different choices of t interval, each of 
these propagators was fitted to the asymptotic form in Eq. ([4.3| ) by minimizing the fit's 
correlated x^- The upper limit of each fitting interval tmax we fixed at the largest t for which 
we had statistically significant propagator data. The lower limit of the fitting interval tmm 
was then progressively increased from 1 to tmax — 2. As t^m was increased, the fitted mass 
and the fit's x^ P^r degree of freedom both generally decreased and the statistical error bar 
increased. For each n and s, the final choice of tmin we took to be the smallest value for 
which the corresponding mass was within the error bars of all the fits with the same n and 
s and larger tmin- Our intent in this procedure was to extract a mass from the largest time 
interval for which the propagator for each combination of n and s was consistent with the 
asymptotic form of Eq. ([4.3|) . 

The solid horizontal lines in Figures - ^ and Figures § - 11^ show the best fitted mass in 
each case and extend over the interval of t on which these fits were made. The dashed lines 
in these figures extend the solid lines to smaller t to show the approach, with increasing t, 
of effective masses to the final mass values. 

For the lattice 16^ x 24 with f3 of 5.93, Tables |TTT| - |VJ show the results of our fits for 
all the combination of n, s, tmin and tmnr which we examined. The best choice of tmin and 



tmax turned out to be 2 and 8, respectively, for all n and s. Tables [VI1| - ^ show the fitted 
masses found with the best choice of tmin and tmax for all the lattice sizes, (3, n and s for 
which our effective mass data showed a plateau at large t. 

As expected, for each lattice size and (3, the fitted masses in Tables |V11| - ^ vary with 
n and s by an amount generally less than the statistical uncertainty in each mass. There is 
also a weak tendency for masses to fall initially with increasing n and s, as the corresponding 
operator's relative coupling to the lightest glueball increases. Then, in some cases, when 
n and s become too large the coupling to the lightest state decreases, the fitted masses 



show some tendency to rise again. To reduce this small remaining statistical uncertainty 
and systematic bias, our final value of mass for each lattice size and (3 was obtained by an 
additional fit of a single common mass to a set of masses from a range of several n and 
s. The common mass was chosen to minimize the correlated x^ of the fit of the common 
mass to the collection of different best mass values. The correlation matrix among the best 
mass values was determined by the bootstrap method. The set of n and s used in each final 
fit was chosen by examining a decreasing sequence of sets, starting with all n and s, and 
progressively eliminating the smallest and largest n and s until a x^ per degree of freedom 
below 2.0 was obtained. The final fit was taken from the largest set of channels yielding a 
X^ below 2.0. If several sets of equal size gave x^ per degree of freedom below 2.0, we chose 



among these the set with smallest x^ per degree of freedom. Tables pCII[ - [XVI| show these 
combined fits and the set of n and s chosen for the final mass value for each lattice and /3. 
In all of these tables, it is clear that once enough of the largest and smallest n and s are 
ehminated to give an acceptable x^ per degree of freedom, the fitted values vary only by 
small fractions of their statistical uncertainty as additional changes are made in the set of 
n and s. The final mass values are collected in Table PCVII. 



At several points in Tables pCII[ - pCVI| , combined fits including several nearby values of n 
and s yield large x^ while separate fits to smaller subsets of n and s give nearly equal masses 
and acceptable x^- This phenomenon, we have found, does not indicate a problem with our 
data or our fits and arises instead because propagators with nearby values of n and s in some 
cases are very highly correlated and yield slightly different masses. A similar problem would 
arise in trying to fit a single value x to, say, a gaussian random variable X with dispersion 
1, and a shifted copy X + 0.0001. For any choice of x the fit's x^ is infinite. Nonetheless, 
for a Monte Carlo ensemble of 1000 X values, taking x as either < X > ±1/a/I000 or 
< X > +0.0001 ± 1/VlOOO is a reliable estimate of the mean of X with systematic error 
much smaller than the statistical error. 

An alternative way to extract a single mass from glueball propagators for a range of n, 
e and s uses the matrix of propagators 

M|'5'-""(ti - 12) =< 5'=^''(ti)5""(t2) > - < S^^'iU) >< ^"^"(t2) > . (4.4) 

For large t and lattice time direction period L, M|''^"^^(t) has the asymptotic form 

Mf ™"(t) ^ Z'=^™"{ea;p(-mt) + exp[-m{L - t)]} (4.5) 

where m is the mass of the lightest scalar glueball and Z^^^"''^^ is a matrix independent of t. 
In principle, M|''''"'"(t) can be extracted from our data and fitted to Eq. ([4.5|) to produce 
a value for m. To find the best m and Z^^^"''^^ by minimizing the fit's x^, however, requires 
the statistical correlation matrix among the fitted M|''™^^(t). If we fit, for example, to three 
choices of {k,S,r), three choices of {n,e,s) and four values of t, the correlation matrix has 
1296 entries. Our underlying data set is too small to provide reliable entries for such a 
large correlation matrix. As a consequence the value of m determined this way will have a 
statistical error which can not be estimated reliably. In practice, we found that the value 
of m produced by this method was not stabile as we varied the sets of {m, 6, r) and {n, e, s) 
and the range of t used in the fit. 



V. TENSOR PROPAGATORS AND MASSES 



A propagator for tensors is defined to be 



Tjnes 



(tl - h 



E 



[< Tnh)Tr{h) >-< rrih) >< Tr{t2) >] 



(5.1) 



where Ti and T2 are the tensor glueball operators of Eq. ( p.4| ), and P^'^*(ti — t2) is then 
averaged over reflections and time direction displacements of ti and t2 f o reduce statistical 
noise. 

Tensor propagators were found for gauge configuration ensembles and operator param- 
eters listed in Tables | and |T[ A tensor glueball mass was extracted from propagators by 
fitting the data to the tensor version of Eq. ( [4.3| ). We obtained a satisfactory tensor glueball 
mass signal only for the lattices with j3 of 5.93, 6.17 and 6.40. We did not find an accept- 
able tensor signal at (3 of 5.70. Overall, the statistical errors in the tensor data are larger 



than those in the scalar data of Section l^ and, as a result, the fitting process encounters 
complications not present in the scalar fits. 

Tables PCVIII| - PCXI] list tensor masses for each gauge ensemble with /3 of 5.93 and above, 
for each set of operator parameters in Table ^ fitted on one or, in some cases, two choices 
of time interval. For all fits the high end of the fitting range tmax is chosen to be the largest 
value at which a statistically significant effective mass is found. The low end of the fitting 
range tmin is then progressively increased. The smallest tmin yielding a mass within one 
standard deviation of the masses for all larger tmin is selected as the lower bound for an 
initial choice of the fitting range. For the lattice 16^ x 24 at /5 of 5.93 and for the lattice 
32^ X 30 X 40 and /3 of 6.40, however, we found that for almost all choices of operator 
parameters a tmin one unit larger than the initial choice yielded a noticeably lower mass. 
These second values of tmin and the corresponding masses are also listed in Tables PCIXJ, and 

ixxn 



Effective mass plots for tensors are shown in Figures 14 - 17, for the four lattices with (3 



of 5.93 and larger, for typical choices of operator parameters. The solid line in each figure 
indicates the mass obtained from a fit over the time interval which the line spans. The 
dashed line in each figures extend the solid line to smaller t to show the approach of effective 
masses to the fitted value. 



Tables PCXIIJ - |XXV| list tensor masses found by combining, as discussed in Section |^, 
the masses fitted to various sets of operators and choices of time interval. Table pCXU 
corresponds to the lattices 12^ x 24 at /3 of 5.93 with fits using the single time interval 



given in Table pCXll| . Table |XX1V| corresponds to the lattice 24^ x 36 at /? of 6.17 with 
fits using the single time interval in Table PCX[ In Tables pCXIl] and PCXIV| , all combined 
fits with acceptable x^ P^r degree of freedom give masses consistent with each other to 
within statistical uncertainties. In each case, the mass corresponding to the largest set with 
acceptable x^, marked with an arrow, is chosen as the final value. 

Table pCXIII| for the lattice 16^ x 24 at /5 of 5.93 shows combined fits using both choices 
of tmin of Table |XIX| . The combined fits using the smaller tmin have unacceptably high x^ 
per degree of freedom. For the fits using the larger tmin the x^ is acceptable, and the fitted 
masses are all consistent with each other within statistical uncertainties. The mass for the 
largest set of operators with the larger tmin is chosen as the final number. Table pCXV| for 



the lattice 32^ x 30 x 40 at /9 of 6.40 also gives combined fits for both tmin in Table pCXl . 
Most fits for both tmin have acceptable x^ P^r degree of freedom. The masses obtained from 
the larger t^m all lie one standard deviation or a bit more below the masses found with the 
smaller tmin, however, and all have significantly better x^ than the fits with the smaller tmin- 
The mass found from the largest set of operators for the larger tmin is therefore chosen as 
the final result. 

The collection of final tensor masses is listed in Table IXXVl. 



VI. VOLUME DEPENDENCE 



We now consider an estimate of the difference between the scalar and tensor glueball 
masses in Table pCVIl] and PCXV1| for finite lattice period L and the infinite volume limits of 
these quantities. 

For large values of L, scalar mQ{L) and tensor m2{L) glueball masses deviate from their 
infinite volume limits, tuq and m2, respectively, by UTT 



ms{L) = ms{l - g, 



exp{- 



- VSniQ L ' 



rrioL 



exp{-moL) 
uioL 



(6.1) 



where s is or 2. In Ref. |T^ for (3 near 6.0, data for mo(-^) is shown to fit the two leading 
terms in Eq. ( |6.1|) reasonably well at 4 values of L ranging from Q/itlq to 12/?7io. This result 
is plausible since for L ranging from G/rriQ to 12/mo, the third term in Eq. ( |6.1| ) is smaller 
than the second by a factor ranging from 0(0.4) to 0(0.2). For our data at (3 of 5.93, 
Table |XVII| shows that mo is above 0.75 so that L of 12 and 16 are larger than 8 /mo and 
12/mo, respectively. Thus we believe that for the data at /3 of 5.93, the leading two terms of 
Eq. ( |6.1|) likely provide a fairly reliable estimate of the L dependence of m,o{L) and m,2{L). 
Fitting the f3 = 5.93 data in Table XVll to the two leading terms of Eq. (|6.1| ) yields mo 
of 0.783 ± 0.012 and go of 1500 ± 1100. In addition a bootstrap calculation yields with 95% 
probability 



mo-mo(16) 



< 0.0037. 



mo 



(6.2) 



. At /? = 5.93, Table pCXVl| combined with the leading two terms of Eq. ( |6.ip gives m,2 of 
1.236 ± 0.037 and g2 is 1300 ± 1200. A bootstrap calculation yields with 95% probability 



m2 - m2(16) 



< 0.0048 



m2 



(6.3) 



Overall, it appears to us safe to conclude that at (3 of 5.93 the difference between scalar 
and tensor masses for L of 16 and their infinite volume limits are of the order of 0.5% or 
less. In Section [VI I| we show that the scalar and tensor glueball masses in Tables pCVII| and 
XXV1| with (3 of 5.93 and greater and m^oL fixed at about 13 are not far from asymptotic 
scaling. We therefore expect the fractional volume dependent errors found in these masses 
to be about the same as the errors at (3 of 5.93. Thus the finite volume errors in all masses 
in Tables pCVII| and [XXVI| with (3 of 5.93 and greater and m,oL of about 13 should be 0.5% 
or less. 



VII. CONTINUUM LIMIT 



The nonzero lattice spacing scalar and tensor glueball masses in lattice units given in 
Tables PCVII| and PCXVI] , respectively, we now convert to physical units and extrapolate to 
zero lattice spacing. 

To convert masses in lattice units to physical units, we divide by a known mass measured 
in lattice units. One natural choice for this conversion factor is the rho mass mp{a)a. Values 



of mp{a)a for three of the four /? in Tables pCVll| and pCXVl| are given in Ref. 0]. For the 



largest /? in Tables PCVII] and |XXVI| , Ref. does not report mp{a)a. For the three /3 consid- 



,(0) 



ered in Ref. [0, however, the ratio [Ajgi^a]/[mp(a)a] is found to be independent of /5 to within 



,(0) 



statistical errors. Here A^g^a is obtained by the 2-loop Callan-Symanzik equation from cnj^g 



found from its mean-field improved [^ relation to /3. Since [Aj^]/[mp{a)a] is constant 



within errors, converting to physical units using A.jrj^a then extrapolating to zero lattice 



(0) 



MS 



a)a. 



Table XXVII 



spacing should give results nearly equivalent to those found using rup 
lists, for each /?, the corresponding mean-field improved aj^ and A^g^a. 

The f3 dependence of valence approximation glueball masses is determined entirely by 
the pure gauge part of the QCD action. The leading irrelevant operator in the pure gauge 
plaquette action has lattice spacing dependence of O(a^). Thus for scalar and tensor glueball 
masses mo and 1712, respectively, we extrapolate to the continuum limit by 



'ms{a)a rUs 



aS^ 



A. 



(0^ 
A/5 



+ C[A 






(7.1) 



where s is or 2. 

If AjTo'^ ill Eq. (|7. 1|) were replaced by mp{a)a, then since the leading irrelevant operator 
in the quark action has lattice spacing dependence of 0{a) it might be argued that the 
quadratic Ola"^) term in the equation's right hand side should be a linear 0{a). This in 
turn would contradict our claim that extrapolation using either mp(a)a or Ai^a will give 
nearly equal results. An answer to this objection is that the approximate constancy of 
[Ajj^a\/[mp(a)a\ implies that the 0{a) irrelevant contribution to mp{a)a is quite small. The 



,(0) 



constancy of [Ajj^]/['mp{a)a\ as a function of a or equivalently as a function of (3 can not 



(0) 



be explained by a cancellation of an 0{a) term in A^^r-a with an 0{a) term in m 



,(0) 



MS 



[a)a 



since A-j-^a is defined to fulfill the true continuum two-loop Callan-Synamzik equation and 
itself has no 0{a) corrections. The leading correction to the /3 dependence of A--j^a is by 



MS 



a multiplicative factor of [1 + 0(/5^)]. If Aj-f^a is replaced by mp{a)a, any significant a 



MS 



dependence which appears will come from the 0(0^) term in ms{a)a. Thus Eq. ( [7.1| ) even 
with mp{a)a substituted for Aj^-a will remain correct. 



(0) 



The scalar data of Tables pCVll| combined with Aj-J—a of Table pCXVll] fitted to Eq. { \l.l\ ) 



at the three largest /5 is shown in Figure |T8[ The predicted continuum limit itiq/AjJ-^ 
is 7.016 ± 0.167. The fit in Figure O has a x^ of 0-6 over a range in which the term 



(0) 



[Aij^a] varies by more than a factor of 3.4. The variation of [m s{a) a] /[Aj^] over the 
fitting range, however, is only slight. Each of the three nonzero lattice spacing values of 
[ms{a)a]/[Aj:^a] is within 1.6 standard deviations of the extrapolated zero lattice spacing 
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result. Thus we believe the extrapolation to zero lattice spacing is quite reliable and would 
expect the predicted continuum mass to be not very different from what would be obtained 
by any other reasonable, smooth extrapolation of the data. 



The tensor data of Tables |XXV1| combined with A^a of Table PCXV11| fitted to Eq. 



MS 



at the three largest /5, the only P for which tensor masses were found, is shown in Figure p 



,(0) 



The predicted continuum limit 1712/ Aj-f^ is 9.65 ± 0.36. The fit in Figure |T9| has a x of 0.8, 



A/5 



while, as before, the term [A^rr^a]^ in Eq. (|7.1| ) varies by more than a factor of 3.4 over the 
fitting range. 

To obtain scalar and tensor glueball masses in units of MeV, we combine the continuum 



,(0) 



,(0) 



limit A^/rUp of 0.305 ± 0.008 0] with rup of 770 MeV to give A^ of 234.9 ± 6.2 MeV. 
The scalar glueball mass becomes 1648 ± 



MeV and the tensor mass becomes 2267 ± 104 



MeV. The continuum limit results are summarized in Table XXVIIl 



For Ajj-/mp we take the value given in Ref. [0] for a lattice with period of about 2.4 
fermi. For the rho mass obtained at /? of 5.7 from a combination of propagators for rho 
operators with smearing parameters 0, 1 and 2, the 2.4 fermi result differs from the result 
for period 3.6 fermi by a bit over one standard deviation. This difference appears to be 
largely a consequence of a slightly poorer separtion of the rho component of the propagator 
from excited state components in the 2.4 fermi rho mass calculation than in the 3.6 fermi 
calculation |2T|]. For the rho operator with smearing parameter 4, which couples more weakly 
to excited states, the difference at (3 of 5.7 between 2.4 fermi and 3.6 fermi predictions is 
much less than one standard deviation. Thus overall it appears to us reasonable to take the 
2.4 fermi calculations as the infinite volume limit, within statistical errors. The continuum 



limit values of Ajj-frUp for the the data combining smearings 0, 
from smearing 4 are nearly identical. 



1 and 2 and for the data 



VIII. COMPARISON WITH OTHER RESULTS 



An independent calculation of the infinite volume, continuum limit of the valence ap- 
proximation to several glueball masses is reported in Ref. [|l^. A second, more recent, 
calculation appears in Ref. 



Tl|| . A comparison of Ref. [^ with the original analysis [^ of 



our results appears in Ref. [2 



The calculation of Ref. [|10| uses the same plaquette action we use but takes a different 



set of glueball operators. The gauge field ensembles of Ref. []T0[ range from 1000 to 3000 
configurations. For the scalar and tensor masses Ref. [IC] reports 1550 ± 50 MeV and 
2270 ± 100 MeV, respectively. The predicted zero lattice spacing masses are not actually 
found by extrapolation to zero lattice spacing, but are obtained instead from calculations at 
f3 of 6.40 of glueball masses in units of the square root of string tension, ^/o^, then converted 
to MeV using an assumed y/a of 440 MeV with zero uncertainty. The uncertainties given in 
the masses are entirely the uncertainties in the (3 of 6.40 calculations of masses in units of 
/o" and are thus missing at least a contribution from the uncertainty in ^/o^. A graph shown 



in Ref. |TD[ suggests that the (3 of 6.40 value of [mo{a)a]/[Ja{a)a] is about 50 MeV below 
the data's zero lattice spacing limit. An additional error of ±50 MeV in the scalar mass 
is therefore proposed in Ref. [|lOl as a consequence of the absence of extrapolation to zero 
lattice spacing. Since [■mo{a)a]/[Ja{a)a] of Ref. JlOl is clearly rising as lattice spacing falls. 
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it does not appear to us that a symmetric error of ±50 MeV an accurate representation of 
the effect of the absence of extrapolation. If the statistical error and extrapolation error in 
the scalar mass are, nonetheless, taken at face value and combined the result is a prediction 
of 1550 ±71 MeV. No estimate is given for the extrapolation error in the tensor mass, which 
is found to be only weakly dependent on lattice spacing if measured in units of ^Jo . A scalar 
mass of 1550 ± 71 MeV is a bit over one standard deviations below the result 1648 ± 58 MeV 
in Table pCXVlll|, while the tensor mass of 2270 ± 100 MeV is in close agreement with our 



value of 2267 ± 104 MeV. 

If the continuum limit of the Ref. |ni| data is found by extrapolation to zero lattice 
spacing of [mo(a)a]/[Ay^a], following Section |V11|, the result for vn^j Isxl^ is 6.67 ± 0.33. 



MS 



Converted to MeV using A^ of 234.9 ± 6.2 MeV, mo becomes 1567 ± 88 MeV. This value 



MS 



are 
found 



is less than a standard deviation below the prediction 1648 ± 58 MeV in Table PCXVIIL 

The calculation of Ref. [O uses an improved action with time direction lattice spacing 
chosen smaller than the space direction. The gauge field ensembles range in size from 
4000 to 10000 configurations. Masses measured in units of the parameter r'^^ [^ 
extrapolated to zero lattice spacing, then converted to MeV using a value of r^^^ 
by extrapolation of r^^ /rrip to zero lattice spacing. As a result of working at relatively 
large values of lattice spacing, some ambiguity is encountered in matching the scalar mass's 
lattice spacing dependence to the small lattice spacing asymptotic behavior expected for 
the improved action. Taking this uncertainty into account, the scalar mass is predicted 
to be 1730 ± 94 MeV. The tensor mass, for which the extrapolation to zero lattice spacing 
encounters no problem, is predicted to be 2400 ±122 MeV. Both numbers are a bit under one 



standard deviation above the predictions in Table pCXVIIl] . For the ratio m2/mQ Ref. |TT 



predicts 1.39 ± 0.04, in good agreement with the value 1.375 ± 0.066 in Table pCXVIII| . Thus 
the difference between Table pCXVlll| and Ref. [^ is almost entirely a discrepancy in overall 
mass scale. 

Combining our extrapolation of 6.67 ±0.33 for the data in Ref. [|r^ with 7.016 ±0.167 in 



Table IXXVlIl gives 6.95±0.15 for mo/A^, thus 1631±55 MeV. Combining 1631±55 MeV 



with 1730 ± 94 MeV of Ref. [^ gives a world average valence approximation scalar mass 
of 1656 ± 47 MeV. This number is consistent with the unmixed scalar mass of 1622 ± 29 
MeV found in Ref. taking the observed states /o(1710), /o(1500) and /o(1400) as the 
mixed versions of the scalar glueball and the two isoscalar spin zero quarkonium states, 
respectively. The state /o(1710) in this calculation is assigned a glueball component of 
73.8 ± 9.5 %. Combining 2270 ± 100 MeV, 2267 ± 104 MeV and 2400 ± 122 MeV gives a 
world average valence approximation tensor mass of 2302 ± 62 MeV. 
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TABLES 



f3 lattice skip count bin 

5^70 le-'^ X 24 50 8,094 8~ 

5.93 12^ X 24 25 48,278 16 

16^ X 24 25 30,640 16 

6.17 243 X 36 25 31,150 16 

6.40 30 X 32^ X 40 25 25,440 16 



TABLE I. Configurations analyzed. 



P lattice n e s 

sTtO 16^ X 24 3-10 025 1-4 

5.93 12^x24 5-7 1.00 4-6 

16^x24 5-8 1.00 3-7 

6.17 24^x36 7,8 1.00 7-10 

6.40 30x32^x40 6,8 1.00 7-11 



TABLE IL Glueball operator parameters. 



n 


s 


*min 


''■max 


mass 


xVd.o.f. 


5 


4 


1 


8 


0.827 ±0.006 


2.26 


5 


4 


2 


8 


0.792 ± 0.013 


1.14 


5 


4 


3 


8 


0.817 ±0.029 


1.22 


5 


4 


4 


8 


0.831 ± 0.068 


1.61 


5 


5 


1 


8 


0.810 ± 0.006 


0.85 


5 


5 


2 


8 


0.791 ±0.013 


0.55 


5 


5 


3 


8 


0.798 ± 0.028 


0.67 


5 


5 


4 


8 


0.767 ±0.059 


0.80 


5 


6 


1 


8 


0.824 ± 0.007 


1.32 


5 


6 


2 


8 


0.793 ± 0.014 


0.58 


5 


6 


3 


8 


0.785 ± 0.031 


0.71 


5 


6 


4 


8 


0.721 ± 0.066 


0.57 



TABLE III. Fitted scalar glueball mass in lattice units for various choices of n, s and fitting 
range for the lattice 16^ x 24 at /3 = 5.93. 
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n 


s 


tmin 


I-maa; 


mass 


xVd.o.f. 


6 


4 


1 


8 


0.815 ± 0.006 


1.82 


6 


4 


2 


8 


0.785 ± 0.012 


0.96 


6 


4 


3 


8 


0.813 ±0.027 


0.90 


6 


4 


4 


8 


0.806 ± 0.064 


1.20 


6 


5 


1 


8 


0.800 ± 0.006 


0.70 


6 


5 


2 


8 


0.785 ± 0.012 


0.50 


6 


5 


3 


8 


0.799 ± 0.027 


0.55 


6 


5 


4 


8 


0.755 ± 0.054 


0.53 


6 


6 


1 


8 


0.811 ±0.006 


0.96 


6 


6 


2 


8 


0.788 ± 0.013 


0.51 


6 


6 


3 


8 


0.789 ± 0.028 


0.64 


6 


6 


4 


8 


0.722 ± 0.060 


0.39 



TABLE IV. Fitted scalar glueball mass in lattice units for various choices of n, s and fitting 
range for the lattice 16^ x 24 at /3 = 5.93. 



n 


s 


tmin 


''■max 


mass 


xVd.o.f. 


7 


4 


1 


8 


0.803 ± 0.006 


1.41 


7 


4 


2 


8 


0.777 ±0.012 


0.70 


7 


4 


3 


8 


0.809 ± 0.026 


0.48 


7 


4 


4 


8 


0.782 ± 0.059 


0.58 


7 


5 


1 


8 


0.792 ± 0.006 


0.61 


7 


5 


2 


8 


0.779 ± 0.012 


0.45 


7 


5 


3 


8 


0.798 ± 0.027 


0.42 


7 


5 


4 


8 


0.746 ± 0.052 


0.25 


7 


6 


1 


8 


0.800 ± 0.006 


0.75 


7 


6 


2 


8 


0.783 ± 0.013 


0.47 


7 


6 


3 


8 


0.793 ± 0.027 


0.55 


7 


6 


4 


8 


0.723 ± 0.056 


0.20 



TABLE V. Fitted scalar glueball mass in lattice units for various choices of n, s and fitting 
range for the lattice 16^ x 24 at /? = 5.93. 
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n 


s 


tmin 


I-maa; 


mass 


xVrf-o./. 


8 


4 


1 


8 


0.789 ± 0.007 


0.97 


8 


4 


2 


8 


0.765 ± 0.013 


0.49 


8 


4 


3 


8 


0.803 ± 0.031 


0.18 


8 


4 


4 


8 


0.754 ± 0.063 


0.04 


8 


5 


1 


8 


0.784 ± 0.006 


0.63 


8 


5 


2 


8 


0.768 ± 0.012 


0.40 


8 


5 


3 


8 


0.791 ± 0.029 


0.31 


8 


5 


4 


8 


0.748 ± 0.055 


0.23 


8 


6 


1 


8 


0.791 ± 0.006 


0.77 


8 


6 


2 


8 


0.775 ± 0.013 


0.57 


8 


6 


3 


8 


0.797 ±0.028 


0.54 


8 


6 


4 


8 


0.716 ± 0.057 


0.04 



TABLE VI. Fitted scalar glueball mass in lattice units for various choices of n, s and fitting 
range for the lattice 16^ x 24 at /3 = 5.93. 



n 


s 


tmin 


''max 


mass 


xVd.o.f. 


4 


1 




5 


0.971 ±0.019 


0.23 


4 


2 




5 


0.959 ±0.017 


0.35 


4 


3 




5 


0.952 ± 0.016 


0.39 


4 


4 




5 


0.984 ± 0.023 


0.26 


5 


1 




5 


0.964 ± 0.018 


0.26 


5 


2 




5 


0.956 ±0.017 


0.29 


5 


3 




5 


0.953 ± 0.017 


0.23 


5 


4 




5 


0.983 ± 0.020 


0.15 


6 


1 




5 


0.958 ±0.017 


0.27 


6 


2 




5 


0.954 ± 0.017 


0.23 


6 


3 




5 


0.956 ±0.017 


0.13 


6 


4 




5 


0.985 ± 0.020 


0.09 



TABLE VIL Fitted scalar glueball mass in lattice units using the best tmin and tmax for various 
choices of n and s for the lattice 16^ x 24 at /3 of 5.70 
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n 


s 


t-min 


trnax 


mass 


xVd.o.f. 


5 


4 


3 


7 


0.752 ± 0.021 


1.30 


5 


5 


3 


7 


0.737 ± 0.020 


1.01 


5 


6 


3 


7 


0.747 ± 0.022 


0.46 


6 


4 


3 


7 


0.754 ± 0.020 


0.97 


6 


5 


3 


7 


0.742 ± 0.020 


0.68 


6 


6 


2 


7 


0.772 ± 0.010 


0.51 


6 


6 


3 


7 


0.751 ± 0.022 


0.29 


7 


4 


3 


7 


0.757 ±0.020 


0.72 


7 


5 


3 


7 


0.747 ± 0.020 


0.42 


7 


6 


2 


7 


0.772 ± 0.010 


0.28 


7 


6 


3 


7 


0.756 ± 0.021 


0.16 



TABLE VIII. Fitted scalar glueball mass in lattice units using the best tmin and tmax for 
various choices of n and s for the lattice 12^ x 24 at /3 of 5.93 



n 


s 


tmin 


''■max 


mass 


xVd.o.f. 


5 


4 


2 


8 


0.792 ± 0.013 


1.14 


5 


5 


2 


8 


0.791 ±0.013 


0.55 


5 


6 


2 


8 


0.793 ± 0.014 


0.58 


6 


4 


2 


8 


0.785 ± 0.012 


0.96 


6 


5 


2 


8 


0.785 ± 0.012 


0.50 


6 


6 


2 


8 


0.788 ± 0.013 


0.51 


7 


4 


2 


8 


0.777 ±0.012 


0.70 


7 


5 


2 


8 


0.779 ± 0.012 


0.45 


7 


6 


2 


8 


0.783 ± 0.013 


0.47 


8 


4 


2 


8 


0.765 ± 0.013 


0.49 


8 


5 


2 


8 


0.768 ± 0.012 


0.40 


8 


6 


2 


8 


0.775 ± 0.013 


0.57 



TABLE IX. Fitted scalar glueball mass in lattice units using the best tmin and tmax for various 
choices of n and s for the lattice 16^ x 24 at /3 of 5.93 
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n 


s 


tmin 


I-maa; 


mass 


xVd.o.f. 


7 


7 


4 


9 


0.570 ±0.018 


0.22 


7 


8 


4 


9 


0.561 ±0.020 


0.14 


7 


9 


4 


9 


0.554 ±0.024 


0.30 


7 


10 


4 


9 


0.540 ± 0.030 


0.38 


8 


7 


4 


9 


0.562 ±0.019 


0.18 


8 


8 


4 


9 


0.551 ±0.018 


0.17 


8 


9 


4 


9 


0.545 ±0.021 


0.10 


8 


10 


4 


9 


0.534 ±0.027 


0.17 



TABLE X. Fitted scalar glueball mass in lattice units using the best tmin and t^ax for various 
choices of n and s for the lattice 24^ x 36 at (3 of 6.17 



n 


s 


tmin 


''■max 


mass 


xVd.o.f. 


6 


7 


4 


12 


0.461 ±0.013 


0.50 


6 


8 


4 


12 


0.446 ±0.012 


0.57 


6 


9 


3 


12 


0.448 ± 0.008 


0.85 


6 


10 


3 


12 


0.435 ± 0.009 


0.57 


6 


11 


3 


12 


0.431 ±0.010 


0.29 


8 


7 


4 


12 


0.456 ±0.013 


0.37 


8 


8 


4 


12 


0.447 ±0.012 


0.54 


8 


9 


4 


12 


0.434 ±0.011 


0.53 


8 


10 


3 


12 


0.433 ± 0.008 


0.61 


8 


11 


4 


12 


0.417 ±0.013 


0.63 



TABLE XL Fitted scalar glueball mass in lattice units using the best tmin and tmax for various 
choices of n and s for the lattice 30 x 32^ x 40 at /? of 6.40 



n 


s 


mass 


xVd.o.f. 


4,5,6 


1,2,3,4 


0.978 ± 0.013 


4.85 


4,5,6 


1,2,3 


0.973 ± 0.014 


3.70 


4,5 


1,2,3,4 


0.975 ± 0.014 


5.39 


5,6 


1,2,3,4 


0.974 ± 0.014 


6.46 


4 


1,2,3,4 


0.965 ± 0.014 


8.95 


5 


1,2,3,4 


0.964 ± 0.015 


9.90 


6 


1,2,3,4 


0.966 ± 0.015 


10.20 


4,5,6 


2,3 


0.973 ± 0.014 


5.07 


5,6 


2,3,4 


0.973 ± 0.014 


7.25 


5 


2,3,4 


0.964 ± 0.015 


14.80 


4,5 


2,3 


0.966 ± 0.014 


3.46 


5,6 


2,3 


0.969 ± 0.014 


6.04 


4,5 


2 


0.957 ±0.015 


1.70 


4,5 


3 


0.953 ±0.015 


0.20 


5,6 


3 


0.953 ±0.015 


1.67 


5,6 


2 


0.955 ±0.015 


0.35 


4 


2,3 


0.955 ± 0.014 


1.36 


5 


2,3 


0.954 ± 0.015 


0.22 


6 


2,3 


0.955 ±0.015 


0.11 ^ 





TABLE XII. Scalar glueball mass in lattice units found by combined fits to sets of n and s for 
the lattice 16^ x 24 at /? of 5.70. The final set chosen is indicated by an arrow. 
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n 


s 


mass 


xVd.o.f. 


5,6,7 


4,5,6 


0.752 ±0.020 


2.35 


5,6,7 


4,5 


0.750 ±0.019 


1.75 


5,6,7 


5,6 


0.747 ± 0.020 


0.98 ^ 


5,6 


4,5,6 


0.752 ±0.019 


3.72 


6,7 


4,5,6 


0.751 ±0.020 


3.67 


5,6 


4,5 


0.749 ±0.019 


2.88 


5,6 


5,6 


0.745 ± 0.020 


1.41 


6,7 


4,5 


0.747 ± 0.020 


2.52 


6,7 


5,6 


0.746 ± 0.020 


1.48 


6 


5,6 


0.742 ± 0.020 


1.19 


7 


5,6 


0.748 ± 0.020 


1.24 





TABLE XIII. Scalar glueball mass in lattice units found by combined fits to sets of n and s for 
the lattice 12^ x 24 at /3 of 5.93. The final set chosen is indicated by an arrow. 



n 


s 


mass 


xVd.o.f. 


5,6,7,8 


4,5,6 


0.781 ±0.011 


1.90 ^ 


5,6,7,8 


4,5 


0.779 ±0.011 


2.77 


5,6,7,8 


5,6 


0.776 ± 0.012 


1.99 


5,6,7 


4,5,6 


0.779 ±0.011 


2.14 


5,6,7 


4,5 


0.778 ± 0.012 


3.24 


5,6,7 


5,6 


0.776 ± 0.012 


2.29 


6,7,8 


4,5,6 


0.782 ±0.011 


2.32 


6,7,8 


4,5 


0.780 ±0.011 


3.50 


6,7,8 


5,6 


0.777 ±0.012 


2.00 


6,7 


4,5,6 


0.778 ±0.011 


3.01 


7,8 


4,5,6 


0.778 ± 0.012 


1.68 


6 


4,5,6 


0.786 ± 0.012 


0.21 


7 


4,5,6 


0.779 ± 0.012 


0.32 


8 


4,5,6 


0.771 ±0.012 


0.86 





TABLE XIV. Scalar glueball mass in lattice units found by combined fits to sets of n and s for 
the lattice 16^ x 24 at /? of 5.93. The final set chosen is indicated by an arrow. 
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n 


s 


mass 


xVd.o.f. 


7,8 


7,8,9,10 


0.559 ±0.017 


1.26 ^ 


7,8 


7,8,9 


0.559 ±0.017 


1.70 


7,8 


8,9,10 


0.553 ±0.017 


1.46 


7 


7,8,9,10 


0.564 ±0.018 


0.58 


7 


7,8,9 


0.566 ±0.018 


0.60 


7 


8,9,10 


0.559 ±0.017 


0.40 


8 


7,8,9,10 


0.555 ±0.017 


0.63 


8 


7,8,9 


0.556 ±0.017 


0.91 


8 


8,9,10 


0.549 ±0.018 


0.30 


7,8 


7,8 


0.559 ±0.017 


2.79 


7 


7,8 


0.566 ±0.018 


1.18 


8 


7,8 


0.556 ±0.017 


1.80 



TABLE XV. Scalar glueball mass in lattice units found by combined fits to sets of n and s for 
the lattice 24^ x 36 at /3 of 6.17. The final set chosen is indicated by an arrow. 



n 


s 


mass 


xVd.o.f. 


6,8 


7,8,9,10,11 


0.4416 ± 0.0074 


3.72 


6,8 


8,9,10,11 


0.4416 ± 0.0074 


4.13 


6,8 


9,10,11 


0.4413 ± 0.0076 


5.37 


6 


9,10,11 


0.4477 ± 0.0082 


9.90 


8 


9,10,11 


0.4363 ± 0.0078 


2.55 


6,8 


10,11 


0.4321 ± 0.0077 


1.14 ^ 


6,8 


10 


0.4327 ± 0.0076 


0.37 


6,8 


11 


0.4281 ± 0.0096 


2.11 


6 


10,11 


0.4360 ± 0.0087 


0.91 


8 


10,11 


0.4326 ± 0.0077 


2.53 



TABLE XVL Scalar glueball mass in lattice units found by combined fits to sets of n and s for 
the lattice 30 x 32^ x 40 at /3 of 6.40. The final set chosen is indicated by an arrow. 
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(5 lattice mass 

5?70 16^ X 24 0.955 ±0.015 

5.93 12^x24 0.747 ±0.020 

16^x24 0.781 ±0.011 

6.17 24^x36 0.559 ±0.017 

6.40 30 X 322 X 4Q 0.4321 ± 0.0077 

TABLE XVII. Final scalar glueball mass values. 



n 


s 


tmin 


I-maa; 


mass 


xVd.o.f. 


5 


4 


2 


5 


1.260 ± 0.036 


1.15 


5 


5 


2 


5 


1.226 ± 0.029 


0.70 


5 


6 


2 


5 


1.224 ±0.031 


0.41 


6 


4 


2 


5 


1.250 ± 0.034 


1.00 


6 


5 


2 


5 


1.218 ±0.029 


0.66 


6 


6 


2 


5 


1.213 ±0.029 


0.42 


7 


4 


2 


5 


1.245 ± 0.032 


0.59 


7 


5 


2 


5 


1.209 ± 0.029 


0.45 


7 


6 


2 


5 


1.206 ± 0.028 


0.31 



TABLE XVIII. Fitted tensor glueball mass in lattice units using the best tmin and tmax for 
various choices of n and s for the lattice 12^ x 24 at /3 = 5.93. 
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n 


s 


tmin 


I- max 


mass 


xVd.o.f. 


5 


4 


1 


5 


1.327 ±0.013 


0.51 


5 


4 


2 


5 


1.284 ± 0.043 


0.34 


5 


5 


1 


5 


1.285 ±0.011 


0.72 


5 


5 


2 


5 


1.270 ± 0.039 


1.02 


5 


6 


1 


5 


1.284 ±0.011 


1.54 


5 


6 


2 


5 


1.252 ± 0.042 


2.04 


6 


4 


1 


5 


1.302 ±0.012 


0.53 


6 


4 


2 


5 


1.282 ± 0.043 


0.70 


6 


5 


1 


5 


1.267 ±0.011 


0.83 


6 


5 


2 


5 


1.264 ± 0.037 


1.24 


6 


6 


1 


5 


1.264 ±0.011 


1.54 


6 


6 


2 


5 


1.243 ± 0.037 


2.18 


7 


4 


1 


5 


1.282 ±0.012 


0.79 


7 


4 


2 


5 


1.278 ± 0.042 


1.19 


7 


5 


1 


5 


1.252 ±0.011 


0.95 


7 


5 


2 


5 


1.255 ± 0.036 


1.42 


7 


6 


1 


5 


1.249 ±0.011 


1.64 


7 


6 


2 


5 


1.233 ± 0.035 


2.37 



TABLE XIX. Fitted tensor glueball mass in lattice units for various choices of n, s, tmin and 
tmax for the lattice 16^ x 24 at /? = 5.93. 



n 


s 


tmin 


tmax 


mass 


xVd.o.f. 


7 


7 


3 


7 


0.861 ±0.029 


0.75 


7 


8 


3 


7 


0.830 ±0.027 


0.24 


7 


9 


3 


7 


0.820 ±0.028 


0.07 


7 


10 


3 


7 


0.813 ±0.034 


0.53 


8 


7 


2 


7 


0.870 ±0.010 


0.46 


8 


8 


3 


7 


0.819 ±0.026 


0.49 


8 


9 


2 


7 


0.839 ± 0.009 


0.27 


8 


10 


3 


7 


0.815 ±0.031 


0.08 



TABLE XX. Fitted tensor glueball mass in lattice units using the best tmin and tmax for various 
choices of n and s for the lattice 24^ x 36 at /3 = 6.17. 



23 



n 


s 


tmin 


I-maa; 


mass 


xVd.o.f. 


6 


6 


4 


9 


0.680 ± 0.043 


0.71 


6 


6 


5 


9 


0.642 ± 0.073 


0.84 


6 


7 


4 


9 


0.660 ± 0.034 


0.65 


6 


7 


5 


9 


0.643 ± 0.063 


0.84 


6 


8 


4 


9 


0.652 ±0.029 


0.33 


6 


8 


5 


9 


0.638 ± 0.059 


0.41 


6 


9 


3 


9 


0.688 ±0.016 


0.58 


6 


9 


4 


9 


0.657 ±0.030 


0.38 


6 


10 


3 


9 


0.675 ±0.016 


0.56 


6 


10 


4 


9 


0.670 ± 0.032 


0.69 


6 


11 


3 


9 


0.660 ±0.017 


0.31 


6 


11 


4 


9 


0.665 ± 0.033 


0.38 


8 


6 


4 


9 


0.658 ± 0.038 


0.46 


8 


6 


5 


9 


0.651 ±0.068 


0.60 


8 


7 


4 


9 


0.636 ± 0.030 


0.50 


8 


8 


4 


9 


0.634 ±0.027 


0.19 


8 


9 


4 


9 


0.637 ±0.026 


0.25 


8 


10 


3 


9 


0.660 ±0.014 


0.35 


8 


10 


4 


9 


0.653 ±0.027 


0.42 


8 


11 


3 


9 


0.646 ±0.015 


0.15 


8 


11 


4 


9 


0.653 ±0.026 


0.17 



TABLE XXI. Fitted tensor glueball mass in lattice units for various choices of n, s, tmin and 
j^ for the lattice 32^ x 30 x 40 at /? = 6.4 



n 


s 


mass 


xVd.o.f. 


5,6,7 


4,5,6 


1.188 ±0.025 


2.13 ^ 


5,6,7 


4,5 


1.207 ±0.028 


1.73 


5,6,7 


5,6 


1.195 ±0.025 


1.33 


5,6 


4,5 


1.206 ±0.029 


1.97 


5,6 


5,6 


1.199 ±0.026 


1.29 


6,7 


4,5 


1.205 ±0.029 


1.90 


6,7 


5,6 


1.207 ±0.027 


0.88 



TABLE XXII. Tensor glueball mass in lattice units found by combined fits to sets of n and s 
for the lattice 12^ x 24 at /3 of 5.93. The final set chosen is indicated by an arrow. 
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n 


s 


train 


i^raax 


mass 


xVd.o.f. 


5,6,7 


4,5 


1 


5 


1.256 ±0.011 


31.87 


5,6 


4,5 


1 


5 


1.260 ±0.011 


44.45 


5,6,7 


4,5,6 


2 


5 


1.234 ± 0.034 


0.60 ^ 


5,6,7 


4,5 


2 


5 


1.252 ± 0.036 


0.48 


5,6,7 


5,6 


2 


5 


1.236 ± 0.036 


0.86 


5,6 


4,5 


2 


5 


1.254 ± 0.036 


0.59 


5,6 


5,6 


2 


5 


1.239 ± 0.035 


1.01 


6,7 


4,5 


2 


5 


1.256 ± 0.037 


0.60 


6,7 


5,6 


2 


5 


1.239 ± 0.034 


1.13 



TABLE XXIII. Tensor glueball mass in lattice units found by combined fits to sets of n and s 
for the lattice 16^ x 24 at /3 of 5.93. The final set chosen is indicated by an arrow. 



n 


s 


mass 


xVd.o.f. 


7,8 


7,8,9,10 


0.856 ±0.010 


2.84 


7,8 


7,8,9 


0.856 ±0.010 


3.83 


7,8 


8,9,10 


0.838 ±0.012 


1.24 ^ 


7,8 


7,8 


0.865 ±0.011 


4.19 


7,8 


8,9 


0.838 ±0.012 


1.72 


7,8 


9,10 


0.840 ±0.012 


0.21 



TABLE XXIV. Tensor glueball mass in lattice units found by combined fits to sets of n and s 
for the lattice 24^ x 36 at /3 = 6.17 The final set chosen is indicated by an arrow. 
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n 


s 


mass 


tmin 


^rnax 


xVd.o.f. 


6,8 


6,7,8,9,10,11 


3,4 


9 


0.655 ±0.012 


1.81 


6,8 


6,7,8,9,10 


3,4 


9 


0.659 ±0.014 


2.04 


6,8 


7,8,9,10,11 


3,4 


9 


0.655 ±0.012 


2.06 


6,8 


6,7,8,9 


3,4 


9 


0.676 ±0.015 


1.52 


6,8 


7,8,9,10 


3,4 


9 


0.659 ±0.013 


2.41 


6,8 


8,9,10,11 


3,4 


9 


0.655 ±0.012 


2.22 


6,8 


6,7,8,9,10,11 


4,5 


9 


0.631 ±0.022 


0.56 ^ 


6,8 


6,7,8,9,10 


4,5 


9 


0.635 ± 0.024 


0.65 


6,8 


7,8,9,10,11 


4,5 


9 


0.632 ±0.021 


0.60 


6,8 


6,7,8,9 


4,5 


9 


0.627 ±0.024 


0.56 


6,8 


7,8,9,10 


4,5 


9 


0.635 ±0.025 


0.73 


6,8 


8,9,10,11 


4,5 


9 


0.630 ±0.022 


0.65 



TABLE XXV. Tensor glueball mass in lattice units found by combined fits to sets of n and s 
for the lattice 32^ x 30 x 40 at /3 = 6.4 The final set chosen is indicated by an arrow. 



/? 



lattice 



mass 



5.93 

6.17 
6.40 



12^ X 24 

16''^ X 24 

24^ X 36 

30 X 32^ X 40 



1.188 ±0.025 
1.234 ± 0.034 
0.838 ± 0.012 
0.631 ± 0.022 



TABLE XXVL Final tensor glueball mass values. 



/3 



a- 



MS 



A^ 



5.700 
5.930 
6.170 
6.400 



0.14557 
0.13180 
0.12183 
0.11407 



0.16612 
0.11444 
0.08265 
0.06177 



TABLE XXVII. For each /3, mean-field improved a-rra and Ai-Z^-o obtained from the 2-loop 



Callan-Synamzik equation. 



MS 
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mo/A^g^ 7.016 ±0.167 



m2/A^ 9.65 ± 0.36 

m2/mo 1.375 ± 0.066 

A^/m„ 0.305 ± 0.008 



mp 770 MeV 

A^ 234.9 ± 6.2 MeV 

MS 

mo 1648 ± 58 MeV 

ma 2267 ± 104 MeV 

TABLE XXVIII. Continuum limit scalar and tensor glueball masses and their conversion to 
MeV. 
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FIG. 1. Scalar effective mass as a function of t for the lattice 16^ x 24 at /3 of 5.93 using the 
smeared operator with n of 5, e of 1.0 and s of 3. 
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FIG. 2. Scalar mass fit and scalar effective mass as a function of t for the lattice 16^ x 24 at /3 
of 5.93 using the smeared operator with n of 5, e of 1.0 and s of 4. 
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FIG. 3. Scalar mass fit and scalar effective mass and fitted mass as a function of t for the lattice 
16^ X 24 at /? of 5.93 using the smeared operator with n of 5, e of 1.0 and s of 5. 
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FIG. 4. Scalar mass fit and scalar effective mass as a function of t for the lattice 16^ x 24 at /3 
of 5.93 using the smeared operator with n of 5, e of 1.0 and s of 6. 
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FIG. 5. Scalar effective mass as a function of t for the lattice 16^ x 24 at /3 of 5.93 using the 
smeared operator with n of 5, e of 1.0 and s of 7. 
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FIG. 6. Scalar mass fit and scalar effective mass as a function of t for the lattice 16^ x 24 at /3 
of 5.93 using the smeared operator with n of 6, e of 1.0 and s of 4. 
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FIG. 7. Scalar mass fit and scalar effective mass as a function of t for the lattice 16^ x 24 at /3 
of 5.93 using the smeared operator with n of 7, e of 1.0 and s of 4. 
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FIG. 8. Scalar mass fit and scalar effective mass as a function of t for the lattice 16^ x 24 at /3 
of 5.93 using the smeared operator with n of 8, e of 1.0 and s of 4. 
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FIG. 9. Scalar mass fit and scalar effective mass as a function of t for the lattice 16^ x 24 at /3 
of 5.70 using the smeared operator with n of 6, e of 0.25 and s of 2. 
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FIG. 10. Scalar mass fit and scalar effective mass as a function of t for the lattice 12^ x 24 at 
P of 5.93 using the smeared operator with n of 6, e of 1.0 and s of 6. 
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FIG. 11. Scalar mass fit and scalar effective mass as a function of t for the lattice 16^ x 24 at 
P of 5.93 using the smeared operator with n of 7, e of 1.0 and s of 5. 
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FIG. 12. Scalar mass fit and scalar effective mass as a function of t for the lattice 24^ x 36 at 
P of 6.17 using the smeared operator with n of 7, e of 1.0 and s of 8. 
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FIG. 13. Scalar mass fit and scalar effective mass as a function of t for the lattice 32^ x 30 x 40 
at P of 6.40 using the smeared operator with n of 8, e of 1.0 and s of 10. 
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FIG. 14. Tensor mass fit and tensor effective mass as a function of t for the lattice 12^ x 24 at 
P of 5.93 using the smeared operator with n of 6, e of 1.0 and s of 6. 
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FIG. 15. Tensor mass fit and tensor effective mass as a function of t for the lattice 16^ x 24 at 
P of 5.93 using the smeared operator with n of 6, e of 1.0 and s of 5. 
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FIG. 16. Tensor mass fit and tensor effective mass as a function of t for the lattice 24^ x 36 at 
P of 6.17 using the smeared operator with n of 7, e of 1.0 and s of 9. 



43 



1.5 



1.0 



0.5 - 





• 








• • 


— 


• • • -r 







• W 1 
1 

• 






- 


- 





0.0 

0.0 2.0 4.0 6.0 8.0 10.0 

t 

FIG. 17. Tensor mass fit and tensor effective mass as a function of t for the lattice 32^ x 30 x 40 
at P of 6.40 using the smeared operator with n of 8, e of 1.0 and s of 10. 
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FIG. 18. The scalar glueball mass in units of A^^-^ extrapolated to zero lattice spacing linearly 
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FIG. 19. The tensor glueball mass in units of A-ij^ extrapolated to zero lattice spacing linearly 
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